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uncertainty of O(Aqcd)- This fundamental difficulty of determining the pole mass to an accuracy 
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bottom quarks can be determined accurately from the O(ma^) quarkonium ground state energy. 
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I. INTRODUCTION 



The pole mass is a central parameter in heavy quark physics but suffers from a renormalon 
caused ambiguity 0, In the heavy quark effective theory (HQET) it is defined formally 
as the heavy quark mass at which the residual mass term vanishes. This definition of the 
pole mass is however incomplete. In general, the parameters and the coefficient functions of 
HQET should be related to the fundamental quantum chromodynamics (QCD) Lagrangian 
through the matching process. In particular, in perturbation theory the pole mass is related 
to the high energy mass such as the MS mass mjj^ (= '"^ms('^ms)) 

oo 

"^pole = "^Ms(l + ' (l) 

n=0 

where = as{mjj^) denotes the strong coupling. The problem with this expansion is that 
it is an asymptotic expansion, with a large order behavior 

(r, -\- ],]] 

Pn = cJ-^^i2(3or [1 + Oil/n)] , (2) 

where Cm is a normalization constant and u = /3i/2/5q, with /3o;/?i denoting the first two 
coefficients of the QCD beta function. Because of the factorial divergence the expansion ^ 
contains an intrinsic uncertainty of order 

"^^'"p(-2^)°^^«^°' 

which means that within the perturbation theory the pole mass cannot be determined to an 
accuracy better than O(Aqcd) even if the MS mass is known. 

On the other hand, one may try to extract the pole mass from the predictions of physical 
observables in HQET. But, as well known, this too does not work since the predictions of 
HQET as well suffer from the renormalon caused uncertainty. As an example, consider the 
inclusive semileptonic B decay rate in HQET, 

r{B X„e.J = ^^Kole)V(«.){l + 0[l/(mSe)1} • (4) 

f{as) contains the QCD corrections to the free b quark decay and can be expanded in 
power series of the strong coupling. This expansion also suffers from the factorially growing 
divergence similar to that in Eq. (j2)) and has an intrinsic uncertainty of 0(AQCD/''^poie) 
[H, 0|- Because of this uncertainty the bottom quark pole mass m^^i^ cannot be extracted to 
an accuracy better than 0{Aqcb), even assuming the decay rate is measured precisely and 
the Cabibbo-Kobayashi-Maskawa matrix element Kib is known. 

Formally, this problem can be resolved by Borel resummation of the divergent series. 
Consider again the pole mass expansion (0). The resummed pole mass can be written as 

mpole = "^BR ± i^m , (5) 

where msR and denote the real and the imaginary parts of the Borel integral, respectively. 
The imaginary part arises from the renormalon singularity that gives rise to the large order 
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behavior Q, and is ambiguous; its sign depends on the integration contour taken. In 
general, an unphysical, ambiguous imaginary part in Borel resummation imphes presence 
of nonperturbative effect and is supposed to be canceled by the latter. In the case of the 
pole mass, however, the ambiguity is spurious in that in physical observables it is canceled 
by a corresponding ambiguity in another divergent series in the coefficient functions or the 
matrix elements 0, For example, in the B decay rate the ambiguous imaginary part 

in the resummed pole mass is canceled by the ambiguity in the resummed /. Thus, as far 
as physical observables are concerned, the introduction of the nonperturbative effect is not 
necessary. 

Since the ambiguous imaginary parts in Borel resummations in HQET are canceled in 
physical observables we can entirely ignore them and keep only the real parts, the 'BR' 
quantities. Under this 'BR' prescription of Borel resummation the pole mass is defined as 

rripoic = mBR • (6) 

Throughout the paper we take this as the definition of the pole mass, unless implied oth- 
erwise. By definition the BR mass, when expanded in the strong coupling constant, has 
exactly the same perturbative coefficients as the pole mass in perturbation theory. 

Although this formal argument resolves the renormalon problem the actual computation 
of the 'BR' quantities is another matter. Since the ambiguity in Borel resummation arises 
from the renormalon singularity, to calculate a BR quantity to an accuracy better than 
O(Aqcd) it is necessary to have a precise description of the Borel transform in the region 
that includes the origin and the renormalon singularity. Within the perturbation theory 
only the first few terms of the power expansion of the Borel transform about the origin are 
known, so describing the Borel transform beyond the immediate neighborhood of the origin 
can be a difficult task, since it would require a precise knowledge on the large order behavior. 

Describing the Borel transform accurately about the renormalon singularity is equivalent 
to having an accurate information on the associated large order behavior. Since the essen- 
tial information in the large order behavior that is missing is the normalization constant 
(residue) its computation will make the computation of BR quantity possible. Fortunately, 
the normalization constant can be calculated in perturbation theory p, |^ , and the residue 
of the pole mass expansion can be computed accurately, within a few percent error, from the 
known next-next-leading order (NNLO) calculations of the expansion 0, 05 0]- This 
accurate computation of the pole mass residue is important, because the residues of other 
divergent series whose leading IR renormalon is cancelled by that of the pole mass can be 
obtained to the same accuracy. 

Once the residue is known the Borel transform can be expanded systematically about 
the renormalon singularity. Then by interpolating this expansion with the perturbative 
expansion about the origin we can obtain an accurate description of the Borel transform in 
the region that contains the origin and the singularity. The interpolation itself is encoded 
in what we call bilocal expansion. 

This idea of bilocal expansion of the Borel transform was successfully applied to the heavy 
quark static potential and the inclusive semileptonic B decay . The Borel resummed 

heavy quark potential, for example, at short and intermediate distances agrees remarkably 
well with lattice calculations. 

So far, the renormalon problem in heavy quark physics was handled essentially by avoid- 
ing it by only considering physical observables. Since physical observables are free from 
renormalons, one can avoid the renormalon problem by replacing the pole mass with a 
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renormalon-free (short distance) mass. Instead, the bilocal expansion allows us to confront 
the renormalon problem directly, and render us infrared sensitive quantities at our disposal. 

In this paper we apply this technique to the perturbative expansion of the heavy quarko- 
nium ground state energy. Our purpose is two fold. First we aim to demonstrate that the 
scheme described above indeed allows us to extract a properly defined pole mass accurately 
from experimental data. It is often claimed that a pole mass cannot be determined to accu- 
racy better than OfA pcn ). We provide here a counterproof. Secondly, we aim to re-examine 
the analysis in Ref. [1^ which extracts from the quarkonium binding energy the top and 
bottom quark masses in the pole mass scheme but without taking into account the renor- 
malon effect. We show that a proper handling of the renormalon results in sizable shifts in 
the extracted quark masses. 



II. THE METHOD 

In this section we give a brief summary of the resummation method of the first infrared 
(IR) renormalon in heavy quark physics using the bilocal expansion. For details we refer 
the reader to ll^ 

Assume a physical quantity A{as) has perturbative expansion 

oo 
n=0 

and a renormalon ambiguity proportional to Aqcd- Then the Borel transform A{b) defined 
by the Borel integral, 

A{as) = -- e-"^^'''A{h)dh, (8) 

has the perturbative expansion about the origin. 



n 



and about the renormalon singularity at 6 = 1/2 of the form 

m = (l + ci(l - 26) + C2(l - 2hf + ■■■) + (analytic part) , (10) 

where C and q are real constants and v same as in 0. The coefficients q can be determined 
by expanding Aqcd in «s, and they depend only on the beta function coefficients. In terms 
of the known four loop beta function the first two coefficients can be found as |J| 

_ /3?-/?o/^2 _ + 4/?o3A/?2 - + m - 2/3f ) - 2/33/^0^ 

Au(3t ' 32z/(z/ - l)/58 • ^ ' 

The "analytic part" in (|1(J|) denotes terms analytic on the disk |6 — 1/2| < 1 about the 
singularity. The radius of convergence of the expansion is expected to be bounded by 
the first IR renormalon at 6 = 1/2 and that of the expansion (jlOj) bounded by the second 
renormalon at h = 3/2. The Borel transforms of the pole mass and the quarkonium binding 
energy are expected to satisfy all these conditions. 
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As noted above, to compute the Borel integral (jHI) to an accuracy better than ©(e-^/^^""'') 
we must have an accurate description of the Borel transform in the region that includes 
both the origin and the first IR renormalon singularity. This can be done with the bilocal 
expansion that interpolates the two expansions Q and (fTUj) . 



A{b) = hm An,m(&), (12) 



where 

^n,m(&) 




M 



l + ^c,(l-26)* 



(13) 



The coefficients hn are to be determined by demanding the bilocal expansion have the same 
perturbative expansion about the origin as 0- The first two coefficients Ci^2 are known, so 
taking M = 2 in (|13|) we then have 

ho = ao-C(l + Ci + C2), 

hi = ai- 2CPo[l - C2 + z/(l + ci + cs)] , 

h2 = a2- 4.CPl[2 + z/(3 + ci - Cs) + //^(l + ci + Cs)] , 

/i3 = as - SCl3l{l + z/)[6 + z/(5 + 2ci - ca) + //^(l + ci + cs)] , etc. (14) 

The interpolating Borel transform A^^yi{h) implements the correct nature of the first renor- 
malon singularity and this allows us to resum the renormalon caused large order behavior 
to all orders. A quick comparison between (fTUj) and (fT^ shows that the expansion about 
the origin in (fT^ . those with the coefficients hn-, simulates the "analytic part" in (fTUj) . 

Clearly, a working bilocal expansion requires computation of the renormalon residue. 
The residue can be calculated perturbatively in a straightforward manner following the 
observation in 0, S|- To compute the residue, notice that 

C = R{\) , (15) 

where 

R{h) = {l- 26)^+"i(6) . (16) 

Now the value of the function R{h) at 6 = | can be evaluated by a series expansion in 
powers of |, since R{h) is expected to be analytic on the disk |6| < 1/2 and is bounded at 
h = 1/2. Even though the evaluation occurs precisely on the boundary of the convergence 
disk, the boundedness of R on the boundary guarantees convergence. Note that for the 
series expansion of R{h) only the Borel transform in perturbative form © is needed. 

The convergence of the power expansion for the residue can be improved by rendering 
the function R to be smoother on the disk \b\ < 1/2. This can be done by pushing away 
from the origin the ultraviolet renormalons and the subleading IR renormalons through an 
employment of a conformal mapping. For the case of the pole mass and the quarkonium 
energy we can use the mapping: 



w = , , (17) 



VTTb+ ^l-2b/3 

which maps the first IR renormalon to w = 1/5 and all others onto the unit circle. In w 
plane the residue can be obtained by expanding R[b{w)] about the origin in power series, 
and evaluating it at w = 1/5. 
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III. RELATION BETWEEN THE POLE MASS AND THE MS MASS 

In this section we show that the relation between the pole mass and the MS mass in the 
case of the bottom and top quark can be given to an accuracy much better than (9(Aqcd)- 
Writing the expansion as 

oo 

rripoie = m^[l + M{as)] , = ^ p„a"+S (18) 

n=0 

the Borel transform Ai{b) of JiA{as) has the first IR renormalon singularity precisely of the 
form (imi) . 

To employ the bilocal expansion for the summation of the series for Ai we first need the 
residue Cm for the pole mass, corresponding to C in Eq. (fTUj) . Following the description in 



oiiowm g tr 

miiiii 



the previous section and using the known first three coefficients 

po = 0.4244, 
Pi = 1.3621 - 0.1055n/, 

I n noi 1 

where n/ denotes the number of light quark flavors, we obtain the residue 



P2 = 6.1404 - 0.8597n/ + 0.0211nJ, (19) 



_ f 0.4244 + 0.1151 - 0.0099 ±0.0080 = 0.5296 ±0.0080, (nj = 5) , . 

™ ~ I 0.4244 ±0.1224 ±0.0101 ±0.0080 = 0.5569 ±0.0080, (n/ = 4) ' ^ ^ 

The errors were obtained by taking the differences between the residues evaluated at NNLO 
and at next-next-next-leading order (NNNLO), with the latter computed using the NNNLO 
coefficients estimated following the method in p^. The good convergence of the residue at 
such a low order perturbation is indeed remarkable, and it can be traced to the renormalon 
dominance in the perturbation expansion of the pole mass. 

Using the values for Cm_above and the coefficients (fT^ we can obtain the interpolating 
Borel transforms ^^0,2(^)5 -^1,2(^)5 and J^2,2{b). With these Borel transforms it is now easy 
to do the actual Borel summation. The integration in the Borel integral can be performed 
easily in w plane, defined by (jl7|) . with the contour now along a ray off the origin to the 
unit circle in the first (or fourth) quadrant 17J. Taking the real part of the integral we have 



' mj^(l ± 0.060336 ± 0.001310 - 0.000029 ± 0.000046) , (n/ = 5 , top quark) 
mj^ll ± 0.1577 ± 0.0041 - 0.0003 ± 0.0001) , {uf = 4 , bottom quark) ' 

(21) 

In this computation we have taken, for demonstration purpose, Cif^{mjj^) = 0.108 and 
a|fl(m(^) = 0.22 for the top and bottom quarks, respectively. The errors are from the 
uncertainties in the computed residues in ()20|). Notice the remarkable convergence compared 
to the unresummed series at the same a^, 

_ f mj^(l ± 0.0458 ± 0.0097 ± 0.0030) , {nj = 5 , top quark) , . 

"^P°'*^ ~ \ mj^(l ± 0.093 ± 0.045 ± 0.032) , (ri/ = 4 , bottom quark) ' ^ ^ 

Taking = 165 GeV for the top and = 4.2 GeV for the bottom quark we see the 
uncertainties in the resummed relations are less than 10 MeV and 1 MeV, respectively — 
much smaller than Aj^. 
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FIG. 1: Interpolating Borel transform Ai2,2{b) (solid) vs ordinary NNLO perturbative Borel 
transform (dashed). 

The drastic improvement in the convergence comes from the proper treatment of the 
renormalon singularity in the Borel transform. For comparison, we plot in Fig. ^ the 
ordinary NNLO Borel transform in power series against the real part of the interpolating 
Borel transform A42,2ib) at Uf = 4. Though they agree on the domain about the origin, as 
should, the deviation becomes obvious at larger values of b, close and beyond the singularity 
at 6 = 1/2. 

The calculation so far clearly shows that direct resummation of the first infrared renor- 
malon is not only feasible but can be efficient. 



IV. RESUMMATION OF THE BINDING ENERGY 

Having established a precise connection between the pole mass and the MS mass we 
proceed to the resummation of the binding energy of the quarkonium IS" state. Since we 
have nothing new to add to the nonperturbative effects or the heavy quark finite width 
effect we shall assume the quarks are stable and the system is completely perturbative. 
Those effects may be added later on to the perturbative result. 

Since the early works by Titard and Yndurain jisl . lol ] there has been many studies on 
the precision calculation of the quarkonium energy within perturbative QCD. The bind- 
ing energy is known to NNLO, and the partial computations of NNNLO energy are also 
incorporated in recent works 11, ^ 2^ 21 . 



The quarkonium IS* energy Mis is given as the sum of the pole mass and the binding 
energy E, 



15 



2m 



pole 



+ E 



2mBR + E- 



BR - 



(23) 
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In perturbation theory the binding energy can be expanded as 

oo 

E = -/i(/i)^£:„(Oa.(/i)"+S (24) 

n=0 

where /z is the renormahzation scale and 

= CFas{fi)mBR, ^=-, (25) 

and Cp = 4/3. The first three coefficients are known 3 and, for mBRtts S> Aqcd, the 
fourth almost 0, [i^, 13' with only the three loop contribution to the static potential 
not calculated. They are given by Sn = ^Pn with 

Po = 1, Pi = 4/3o^ + h , 

Pa = l2p^L^ + (-8/?o' + 4/?i + 6Poki)L + k2 , 

Pa = 32(3^ + (-56/?o' + 28/5o/?i + 24/52^1)^2 

+(16/?o3 - i6/?o/3i + 4/?2 - I2f3^h + Qf^ih + 8/?oA;2)i^ + h , (26) 
where L = In(^) and 

= 1(97/6 -lln./g), 

TT 

= —(337.9471 -40.9648571^ + 1. 162857^2.), 

1 



^3 = —(7078.7900 - 1215.5475n/ + 69.45081671^ - 1.21475n} 



+0.031250 as + 474.2893 ln[a,(/i)]) , (27) 
and the beta function coefficients, 

/?o = ^(11 - 2n^/3), A = ^(102 - 38n;/3) , 

/?2 = --l-r(2857/2- 5033r2f/18 + 325^2/54). (28) 
(47r)'^ ■' 

03 in ()27p denotes the unknown three loop coefficient in the MS scheme of the static potential 
in momentum space. Using the Borel transform method of estimating unknown higher order 
coefficients fli^ , we can find an estimate of the three loop coefficient V3 of the static potential 
in coordinates space, 

V{r) = -J2Vna,{l/rr+\ (29) 



^ n=0 



which reads 



r -19.33 + 2.73, inj = 5) 
^3 ~ 1 -27.03 + 3.50 , {nj = 4) " ^"^^^ 



The corresponding 03 is then given by ^ 

03 _/ 34 + 63, {uf = 5) 



43 1 59 + 81 , {uf = 4) ■ '^^^^ 



^ This agrees with the existing estimates within the errors, = 60,98 from the Pade Esj and 

37, 72 from the large order behavior 0, for Uf — 5, 4, respectively. 
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Note the errors in V3 estimates get amplified in the conversion from V3 to 03 by an order 
of magnitude, resulting in large errors in the latter. So a few 100% error in 03 should be 
counted as normal. 

The binding energy has a renormalon caused ambiguity proportional to Aqcd which is to 
be canceled by that of the pole mass term (2mpoie) jl?, IE 23; Borel transform 



E of the binding energy can be expanded about the renormalon singularity as 



:i -26)n 



(1 + 



ci(l - 26) + C2(l - 2bf + ■■■) + (analytic part) (32) 



with Cm and ci,2 given by and (fTTj) . respectively. Note that in Eq. (jH^ the renor- 
malon residue of E is expressed in terms of the pole mass residue Cm, using the renormalon 
cancellation between the pole mass and the binding energy. This is to utilize the accurate 
computation of the pole mass residue. One can, of course, proceed without refereeing to it, 
by directly computing the residue from the coefficients ^I^ : however, the convergence is not 
as good as in the pole mass residue [§, Hi- Combining (jH^ with the expansion about 
the origin, 

n=0 



Jo, 



(33) 



into a bilocal expansion following the description in Sec. |n]we can obtain, using the known 
coefficients fOUj) . the first four interpolating Borel transforms E]y^2{b,^) {N = 0, 1,2,3). 

The Borel resummed binding energy E-br can then be obtained by substituting E{b,C,) 
with i?Ar 2(6, the integral 



E- 



BR 



Re 



2 /"OOztie 



e-''/l3oMf^)E{b,^)db 



(34) 



In the following we take af^{Mz) = 0.1172 ± 0.002 and use the four loop beta 
function for the computation of the running coupling ^^(/i). 



A. Top quark 

Because of the large mass and decay width of the top quark the toponium IS* state is 
expected to have small nonperturbative effects and to be well described by perturbation 
theory. In future e'^e~ linear colliders the toponium energy Mis is expected to be measured 
precisely from the peak position of top threshold production cross section. This provides a 
unique opportunity for precision determination of the top mass. 

We now perform the Borel integration in ()34|) with the interpolating Borel transforms. 
The results are plotted in Fig. |21at the pole mass value ttibr = 175 GeV. 

Notice the very small renormalization scale dependence in the NNNLO resummed energy. 
It is less than 10 MeV over the range 30 GeV < fi < 150 GeV with a^/i^ = 34. For 
comparison we plot the unresummed binding energies in Fig. lousing the same pole mass. 
Compared to the unresummed binding energies the improvement of the resummed is clear. 
They have much smaller perturbation order dependence and scale dependence. It is evident 
that for precise determination of the binding energy the renormalon must be taken into 
account properly. 
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FIG. 2: The resummed toponium 15 state binding energy using i?i^2 (dotted), £^2,2 (dashed), and 
-^3,2 (soUd). The shaded band denotes the variation due to the uncertainty in the estimate of 03 
in Eq. (|5T|) . The central hne, upper, and lower boundaries correspond to 03/4^ = 34, —29, and 97, 
respectively. 
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FIG. 3: The unresummed toponium 15 state binding energy at NLO(dotted), NNLO(dashed), 
and NNNLO (solid). For comparison the resummed binding energies using £^2,2 (short-dashed) and 
-E'3,2 (dot-dashed) are given. The same pole mass was employed in the resummed and unresummed 
cases. 



The dominant theoretical uncertainty in the resummed binding energy comes from the 
uncertainty in the strong coupling constant. In the pole mass scheme under consideration 
there is a natural scale for the binding energy, the Bohr scale /i(/iBohr) = A^Bohr, with /^Bohr ~ 
33 GeV. We expect an optimal result would come around this scale. Thus, taking the 
variation of the binding energy at the Bohr scale under 5af^{Mz) = 0.002 we estimate the 
error caused by the strong coupling to be about ±150 MeV to the binding energy, which 
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FIG. 4: The resummed T(15) binding energy using £Ji ^2 (dotted), £^2,2 (dashed), and £^3^2 (sohd). 
The shaded band denotes the variation due to the uncertainty in the estimate of 03 in Eq. ()31() . 
The central hne, upper and lower boundaries correspond to 03/4^ = 59, —22, and 140, respectively. 

would give ±75 MeV uncertainty in the pole mass. Other than the strong coupling the main 
source of uncertainty is 03, as can be seen in Fig. |21 Comparatively, the renormalization 
scale dependence is small. In Fig. |21we see that the binding energies at next-leading order 
(NLO), NNLO and NNNLO come close together in the range /iBohr < < 80 GeV, and also 
the principle of minimal sensitivity (PMS) [j^H scales for the NNLO and NNNLO curves 
lie within the same interval. The renormalization scale dependence over this interval of 
the central NNNLO line is negligible, less than 2 MeV. The uncertainty in the renormalon 
residue at = 5 in (^01) causes ±5 MeV error to the binding energy. The uncertainty in 03 
causes another ±16 MeV error in the binding energy. Combining these errors we conclude 
that the current theoretical error, excluding that caused by the strong coupling constant, 
on the pole mass is ±10 MeV. 

Once the pole mass is extracted from the binding energy the top mass in the MS scheme 
can be determined using the relation between the BR mass and MS mass discussed in Sec. 
mil The conversion introduces another error due to the uncertainty in the relation, about ±7 
MeV to the MS top mass. Combining these uncertainties in quadrature the corresponding 
theoretical uncertainty in the MS mass is about ±12 MeV. 

B. bottom quark 

The resummation of the T(15') binding energy can proceed in the same way as in the 
toponium, only with a few changes of the parameters. We take the renormalon residue and 
03 estimates in 1^ and (jHTj) . respectively, for = 4, and for the strong coupling af\jj) 
we use RunDec |3^ with the four loop beta function and three loop matching. 

The resummed binding energies from the interpolating Borel transforms at NLO, NNLO, 
and NNNLO are plotted in Fig. E] at the pole mass = 4.9 GeV. For comparison 

the unresummed binding energies are plotted in Fig. |3| Again, the improvement of the 
resummed energy is clear. It has better convergence and less scale dependence. 
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FIG. 5: The unresummed T{IS) binding energy at NLO(dotted), NNLO(dashed), and NNNLO 
(solid). For comparison the resummed binding energies using £'2,2 (short-dashed) and £'3^2 (dot- 
dashed) are given. The same pole mass was employed in the resummed and unresummed cases. 

Before the BR mass can be extracted from the fitting of the resummed energy to the 
precisely measured T(IS') mass, the nonperturbative effects and the charm mass effect should 
be subtracted from the T(IS') mass. Without taking these effects into account we obtain 
mBR ~ 4.9 GeV. The leading nonperturbative effect comes from the Stark effect of the long 
range gluon fields in the QCD vacuum [^HHl- When the correlation time of the quarkonium 
state is much smaller than 1/Aqcd the nonperturbative effect may be expressed in terms of 
local gluon condensate. Whether T(IS') is such a state is not clear, but assuming so, the 
nonperturbative effect is given by 



The gluon condensate is poorly known. In general, fitted power corrections are contaminated 
by perturbative corrections, and sensitive to the perturbation order employed in fitting. For 
gluon condensate, it has been observed [s^ that with higher order perturbation the value 
extracted is much smaller than the original sum rule value in [s^ . Since in the binding energy 
the perturbative correction is already at relatively high order (NNNLO) ^ the condensate 
for Eq. should not contain significant perturbative corrections, and thus likely to be 

smaller than the original sum rule value. Here we take < asG^j, >= 0.02 ± 0.02. Since 
the local condensate approximation for the nonperturbative effect is likely to provide an 
upper bound for any potential nonlocal effect [38] , it may be reasonable to assume that this 
value covers adequately the nonperturbative effect. With mBR = 4.9 GeV we then get the 
nonperturbative effect 60 ±60 MeV in the binding energy. For the nonzero charm mass effect 
we find, using the result in [39.. ,40] , —20 ± 15 MeV in the binding energy. Subtracting these 
two effects from the experimentally measured T(IS') mass and comparing the subtracted 



^ The perturbative ultrasoft correction itself starts at 0{ma^g) [3] 



1872mBR7r < asGl^ 
1275 {Cpasm^^Y 



> 



(35) 
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mass with the PMS value of the resummed energy, 2mBR + i?BR, we obtain meR = 4.89±0.03 
GeV. 

For the error estimate we notice in Fig. EJthe resummed energies at NLO, NNLO, and 
NNNLO converge around /i ^ 2.5 GeV and the PMS scales for the NNLO and NNNLO 
resummed energies lie within the range meohr < < 5 GeV, where mBohr ~ 2 GeV. We 
therefore take the variations of the resummed energies in this range as the error estimate. 
Then the uncertainties in the and the strong coupling give ±27 MeV and ±13 MeV, 
respectively, to the binding energy, and the renormalization scale dependence and the un- 
certainty in the residue (j^UI) cause another ±30 MeV and ±12 MeV, respectively. The error 
due to the unknown higher order contribution should be partly covered by the error from 
the residue. However, setting this fact aside we shall assign an independent error ±35 MeV 
to the binding energy, which is the difference between the respective energies at NNLO and 
at NNNLO at the Bohr scale. It is interesting to observe that the PMS values at NNLO and 
NNNLO are virtually identical, an indication that the unknown higher order effect should 
be small. 

Combining these errors (including those from the gluon condensate and the nonzero charm 
mass effect) in quadrature we have the BR mass for the bottom quark 

m^Si = 4-89 ± 0.04 GeV . (36) 



Now, from the BR mass the MS mass can be obtained through the relation, m^{l + 
7Mbr[q;|^'('^^)]} — ""^BR) from which we obtain 



MS 

Compare this to the result from Ref. 1^ |. 



mg^ = 4.20 ± 0.04 GeV . (37) 



mi^ = 4.346 ± 0.070 GeV , (38) 

which was obtained without taking the renormalon into account. Notice that the central 
values have a significant difference, and our value has a smaller error. The reduction in 
the error largely comes from the smaller renormalization scale dependence. The assigned 
error for the scale dependence in (j38p is more than four times larger than ours, because of 
much stronger scale dependence in its derivation. This comparison shows that renormalon 
effect is significant and must be taken into account for accurate determination of the bottom 
quark mass. Our value for the bottom quark mass ()37|1 compares favorably with the values 



collected in 41 



V. DISCUSSION 

In this paper we have shown that an accurate calculation of the normalization constant 
of the large order behavior and information on the renormalon singularity allows us a pre- 
cision calculation of the pole mass and quarkonium binding energy. The existing method of 
handling the renormalon problem was bypassing it, by concentrating exclusively on directly 
measurable observables and employing renormalon-free short distance masses. In contrast, 
in our approach each of the first IR renormalon in IR sensitive quantities is resummed inde- 
pendently, and as a consequence those IR sensitive quantities can be extracted accurately 
from experimental data. 
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One of the advantages of our approach is that it provides a natural solution to the scale 
mixing problem of the approaches based on renormalon cancellation. In systems with far- 
separated multiple scales the renormalon cancellation approach inevitably mixes those scales 
through the renormalization scale, forcing one to choose an optimal scale over a wide range 
of scales. Depending on the problem this can be a source of significant error. For instance, 
in the top-pair threshold production the normalization of the production cross section is 
known to be sensitive to the choice of the renormalization scale 42], which can limit the 
accuracy of the extracted strong coupling constant. With our approach, this problem can 
be solved by resumming the top pole mass and the interquark potential at their respective 
optimal scales, namely, the top mass and the interquark distance. 

This suggests a potential improvement in the bottom quark mass determination. Instead 
of applying the resummation to the perturbatively calculated binding energy, as we did in 
this paper, one may solve the Schroedinger equation, which was advocated in |43|], with a 
resummed interquark potential to obtain a more accurate binding energy. It should be noted, 
however, presently the error in the bottom quark mass is dominated by the nonperturbative 
effect. 
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